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摘 要 : 开展 干旱 预测 是 有 效应 对 干旱 风险 的 前 提 基 础 ,根据 1960—2016 年 三 江平 原 7 个 站 点 逐 
日 降水 和 气温 数据 ,利用 ARIMA 和 ANN 模型 对 不 同时 间 尺 度 标准 化 降水 燕 散 指数 (SPE1) 序 列 进 
行 分 析 建 模 预测 。 借 助 相关 系数 尺 \、 纳 什 效 率 系数 NSE, Kendall 秩 相关 系数 7、 均 方 误差 MSE 和 
Kolmogorov-Smirnov (K-S) 检验 对 模型 的 有 效 性 进行 了 判定 ,然后 分 别 用 ARIMA 和 ANN 模型 进行 
12 步 预 测 , 并 将 预测 值 与 实际 值 进行 比较 。 结 果 表 明 :(1) ARIMA 模型 和 ANN 模型 对 SPEI 的 预 
测 能 力 都 随时 间 尺 度 的 增加 而 逐渐 提高 。(2) 两 种 模型 对 3.6 个 月 尺度 SPEI 的 预测 精度 偏 低 ,9、 
12 .24 个 月 的 SPEI 的 预测 精度 在 70% 以 上 ;(3) SPEI -9,SPEI - 12,SPEI -24 三 个 时 间 尺 度 ANN 


模型 的 预测 精度 优 于 ARIMA 模型 。 
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干旱 是 人 类 面临 的 主要 自然 灾害 之 一 ,而 我 国 
几乎 每 年 都 会 发 生 ,平均 2 ~3 a 就 会 发 生 一 次 严重 
的 干旱 灾害 ”1 。 加 强 干 旱 监测 .预测 方面 的 研究 
对 相关 部 门 预防 干旱 灾害 ,减少 干旱 灾害 损失 具有 
重要 的 意义 。 目 前 用 以 预测 干旱 的 方法 有 很 
多 ,如 人 工 神经 网 络 (ANN ) 模型 .差分 自 回归 移动 
平均 (ARIMA ) 模型 .支持 向 量 机 (SVM ) 等 。 其 中 ， 
ARIMA 模型 作为 传统 的 时 间 序 列 模型 ,应 用 最 普 
遍 , 应 用 范围 比较 广泛 ,如 韩 萍 等 通过 对 关中 地 
区 多 尺度 标准 化 降水 指数 的 ARIMA 建 模 ,指出 该 
模型 在 干旱 预测 方面 具有 良好 的 效果 ,ZHANG Y 和 
LI W! 利用 ARIMA 等 模型 对 中 国 海河 北 系 的 干旱 
情况 进行 了 预测 分 析 , 说 明了 ARIMA 模型 具有 较 
为 优良 的 预测 精度 。ANN 模型 作为 处 理 非 线性 问 
题 的 典型 模型 ,受到 了 众多 学 者 的 青睐 ,也 越 来 越 多 
的 被 运用 到 各 领域 的 预测 中 ~"  。 其 中 在 干旱 预 
测 领 域 ,BARUA 等 “为 提高 水 资源 管理 效率 ,利用 
ANN 神经 网 络 对 干旱 进行 了 预测 , BELAYNEH 
等 ( 引 利 用 ANN 模型 在 阿 瓦 什 河 流域 对 SPI 指数 进 
行 了 分 析 建 模 , 二 者 都 证 明了 ANN 模型 预测 精度 的 
DER TE, 8] 36 7 JR] FI ANN 人 工 神经 网 络 对 淮河 


© KAH: 2018-05-02; 修订 日 期 : 2018 -07-11 


流域 土壤 水 分 进行 反 演 研究 ,指出 该 模型 对 淮河 流 
域 的 干旱 监测 具有 良好 的 使 用 价值 与 应 用 前 景 。 
干旱 程度 常 利 用 干旱 指数 进行 定量 评价 51 。 
国内 外 常用 的 分 析 指 数 有 : Palmer 干旱 指数 (PD- 
SI) .标准 化 降水 指数 (SP1) .综合 气 象 干旱 指数 
(CI) Aie fe EK Z BIG AC (SPEI) SE, SPEI 融合 
了 SPI 计算 简单 及 多 时 间 尺 度 的 特性 和 PDSI X128 
散 需 求 变 化 的 灵敏 性 等 方面 的 优点 ,是 表征 干旱 程 
度 的 理想 工具 '“ "1。SPEI 具有 多 个 时 间 尺 度 属 
性 ,如 3 个 月 时 间 尺 度 的 SPEI 和 季节 降水 估计 相关 
性 较 高 ;12 个 月 时 间 尺 度 的 SPEI 对 于 下 层 土壤 水 
分 和 河流 径流 量 等 有 较 好 的 反映 ;24 个 月 时 间 尺 度 
B SPEI 可 以 更 深入 地 涉及 到 地 下 水 位 水 库 蓄 水 量 
及 民生 问题 . SPEI 多 时 间 尺 度 的 特征 使 其 在 干 
旱 分 析 、 评 估 等 领域 被 广泛 使 用 ”|。 
目前 的 文献 大 多 以 年 或 某 一 月 份 时 间 尺 度 的 
SPEI 为 研究 对 象 , 而 关于 模型 对 不 同时 间 尺 度 
SPEI 的 适用 性 研究 却 鲜 见 。 鉴 于 此 ,本 文 应 用 三 江 
平原 1960 一 2016 年 7 个 气象 台 站 的 常规 气象 观测 
资料 ,以 不 同时 间 尺 度 (3 个 月 .6 4H 2 个 月 、24 
个 月 ) 的 SPEI 为 基础 ,研究 比较 ARIMA 和 ANN 模 
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型 在 干旱 预测 中 的 适用 性 ,并 对 预测 能 力 进 行 了 评 
价 分 析 ,以 期 能 为 三 江平 原 的 抗旱 减灾 工作 提供 一 
定 的 依据 。 


1 研究 方法 


1.1 ARIMA 模型 

ARIMA 模型 于 1976 年 由 BOX 等 提出 ,是 处 理 
时 间 序 列 的 最 普遍 、 最 直接 的 建 模 方法 ,也 是 对 时 间 
序列 进行 预测 的 科学 方法 。 它 将 预测 对 象 随时 间 变 
化 形成 的 序列 看 作 是 一 个 随机 序列 ,即时 间 序 列 依 
赖 于 时 间 + 的 一 组 随机 变量 。 其 中 ,单个 序列 值 出 
现 具有 不 确定 性 ,但 整个 序列 的 变化 则 呈现 出 一 定 
的 规律 性 ,可 以 用 相应 的 数学 模型 加 以 近似 描述 ,这 
就 是 ARMA 模型 的 基本 思想 ARMA 模型 有 三 
种 类 型 : 自 回 归 ( AR ) 模 型 .移动 平均 (MA ) 模型 以 
及 自 回归 移动 平均 (ARMA ) 模 型 。 只 有 平稳 的 时 间 
序列 才能 够 直接 建立 ARMA 模型 ,现实 中 常见 的 时 
间 序 列 多 具有 某 种 趋势 ,但 往往 通过 差分 可 使 其 满 
足 平稳 性 要 求 。 然 后 则 可 以 建立 ARMA 模型 , 即 
ARIMA 模型 。 通 常 , 经 过 d 次 差分 后 的 ARMA 模型 
方程 由 下 式 表示 : 


Y, =PI] +Y, 2 + wp Visp te, -0,e,., 2m 


05e, -2 eene, (1) 


WE ARIMA(p,d,q) 。 
FINEER TET. B, 则 ARIMA (p,d,q) Zr fé n] 
简 记 为 


9,(B)V'Y,=0,(B)e, (2) 

其 中 ， 
p,(B) =1 -ọ,B -p,B° -9,B" (3) 
0,(B) 1-0,B - 6B. - 6, B' (4) 


Y, JUR 90 [8 p: (1 2 1,2, p) R16, (j 21,2, 
…,g) 分 别 为 自 回 归 系 数 和 滑动 平均 系数 ;e, Du E 
声 序列 , 且 e, - N(0,07) o 

本 文 将 1960 一 2005 年 的 数据 作为 训练 数据 ,由 
数据 序列 的 自 相关 图 和 偏 自 相关 图 确定 阶 数 P 和 4 
的 取 值 范围 ,然后 根据 AIC 准则 对 p RI q 进行 优选 ， 
即 在 不 同 p 和 9g 组 合 的 模型 中 ,选择 最 小 的 AIC 值 
所 对 应 的 参数 ,从 而 得 到 最 优 ARIMA 模型 。 然 后 
利用 2006—2015 年 的 数据 对 选 定 的 最 优 模型 进行 
预测 验证 ,最 后 向 前 预测 了 2016 年 的 SPEI 值 并 将 
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其 与 实际 值 进行 比较 。 这 一 过 程 在 R 语言 平台 实 
现 。 
1.2 ANN 模型 

ANN 人 工 神经 网 络 是 一 种 模仿 动物 神经 网 络 
行为 特征 进行 分 布 式 并 行 信 息 处 理 的 算法 , 按 是 否 
含有 延迟 或 反馈 环节 可 以 分 为 静态 神经 网 络 和 动态 
神经 网 络 。 静 态 神 经 网 络 是 无 反馈 无 记忆 的 ,输出 
仅 依 赖 于 当前 的 输入 ,例如 BP 神经 网 络 和 RBF 神 
经 网 络 ;动态 神经 网 络 ” 是 带 有 反馈 的 ,其 输出 不 
仅 依 赖 当前 的 输入 ,而 且 依 赖 之 前 的 输入 ,通过 反 
馈 , 神 经 网 络 能 将 前 一 时 刻 的 数据 保留 ,使 其 加 入 到 
下 一 时 刻 数据 的 计算 ,使 网 络 不 仅 具 有 动态 性 ,而 且 
保留 的 系统 信息 也 更 加 完整 。 动 态 神经 网 络 具 有 很 
强 的 学 习 能 力 和 通 近 任意 函数 的 特点 ,在 计算 效率 
和 预测 精度 等 方面 都 优 于 静态 神经 网 络 。 自 20 t 
纪 80 年 代 末 以 来 ,动态 神经 网 络 作为 一 种 轩 新 的 非 
线性 模型 被 广泛 地 用 于 变形 预报 .模式 识别 .最 优化 
决策 等 方面 “1!。 因 此 ,本 文选 择 动态 神经 网 络 
(NAR 回归 神经 网 络 ) 进 行 SPEI 时 间 序 列 预测 。 

NAR 神经 网 络 模型 可 以 表示 为 : 

y(t) Sfiy(t - 1) ,y(t -2), ,y(t-d)) (5) 
式 中 :d 为 延 时 阶 数 ; /(， ) 为 神经 网 络 模型 ,可 以 
看 出 y(i) 值 是 由 y(t 71) ,y (£72) y GO -) BS 
值 所 决定 的 。 表 明 该 模型 用 历史 的 数据 值 来 推断 当 
前 及 未 来 的 数据 值 ,具有 动态 性 和 延展 性 。 

NAR 动态 神经 网 络 模 型 的 训练 方法 通常 使 用 
Levenberg-Marquardt 算法 , 因 其 训练 速度 极 快 ,对 于 
中 小 型 网 络 而 言 ,此 训练 算法 是 最 好 的 算法 。LM 
算法 的 计算 公式 为 : 

X,,, 2 X, - (JJ +u) J"e (6) 
WP X, Xi ,1 分 别 是 第 和 hk+1 次 迭代 的 网 络 权 值 ; 
为 雅 可 比 矩 阵 ; 了 为 单位 矩阵 ;uw 为 神经 网 络 学 习 
率 ;e 为 随机 误差 。 

本 部 分 将 1960 一 2015 年 数据 划分 为 三 类 :训练 
数据 、 验 证 数据 和 测试 数据 , 比例 设置 为 :70% 、 
1596 15% ,之 后 设置 神经 元 个 数 、 延 迟 阶 数 等 参数 ， 
并 通过 调节 参数 判断 误差 ,使 得 Autocorrelation 在 满 
意 值 范围 内 ,筛选 出 相对 最 优 的 网 络 模型 ,然后 利用 
训练 成 功 后 的 模型 ,预测 2016 年 的 12 个 值 。 这 一 
过 程 通过 Matlab 软件 实现 。 

1.3 标准 化 降水 蒸 散 指数 (SPET) 
采用 具有 多 时 间 尺 度 特征 的 SPEI 作为 气候 干 


杨 慧 荣 等 ;ARIMA 和 ANN 模型 的 干旱 预测 适用 性 研究 


湿 状 况 的 表征 。 该 指数 计算 原理 简单 易 行 , 即 通过 
衡量 降水 量 和 湾 在 葵 散 量 差 值 偏离 平均 状态 的 程度 
来 判定 干 湿 状况 所” 。 具 体 计 算 步骤 如 下 : 

(1) YFTETE TE Z& BE (PET) , Z& SCÓK H B8 2: 
Thormthwaite 77 12: ?*. 。 


A 
PET - eK x (Ar) (7) 


式 中 :K 为 根据 纬度 计算 的 修正 系数 ;7 为 月 平均 温 
度 ; 刀 为 年 热量 指数 ;4 为 以 五 为 基础 的 系数 。 

(2) 计算 逐 月 降水 与 藻 散 的 差 值 : 

D, =P, - PET, (8) 

rp LP; 为 逐 月 降水 量 ; PET, 29 H ETE ZR CURE E 
(mm), 

(3) 采用 3 个 参数 的 log-logistic 概率 分 布 拟 合 
所 构建 的 数据 序列 D; : 


B. -1 
- a 
F(x) [1 ] (9) 
其 中 ,参数 B voy 采用 工 一 和 矩 估 计 方 法 获得 : 
B (wo -2o,)B 
“ETA +1/B8)T( 1/8) (10) 
2 17 o 
de eris. an 


y 20s - al (C1 * 178) T (1 - 1/8) (12) 


其 中 ,为 阶乘 函数 ,oo o, o 为 数据 序列 D; 的 概 
KIMIE ,计算 方法 如 下 : 


Te (13) 


式 中 :NN 为 参与 计算 的 月 份 个 数 。 
(4) 对 累积 概率 密度 进行 标准 化 ,获取 相应 的 
SPEI fË: 


P-1-F(x) (14) 
当 P<0.5 时 ， 
W - / -2lnP (15) 
Wte W? 
SPEI = W - — L E te (16) 


1+d W+d, W +d, W 


式 中 :cv =2. 515 517 ,c, =0. 802 853 ‚c, =0. 010 328; 
d, =1. 432 788 ,d, =0. 189 269,d, =0. 001 308, 24 
P >0.5 时 , 本 文通 过 划分 SPEI 的 等 级 来 表征 干旱 
情况 ,如 表 1。 


R1 SPEI 干旱 等 级 的 标准 划分 表 
Tab.1 Drought classifcation of SPEI 


SPEI 数值 早 等 级 
-0.5 <SPET<0.5 无 旱 
-1.0 «SPEI« -0.5 轻微 干旱 
-1.5 < SPEIS -1.0 中 度 干旱 
-2.0 < SPEIS -1.5 严重 干旱 
SPEI< -2.0 极端 干旱 
W= /-2ln(1- P) (17) 
c, t c, W E c, W 
SPEI = o 2 w — (18) 


1+d, W+d, W +da W 


1.4 评价 验证 指标 
1.4.1 相关 系数 (R) 


Xi -7) (J; =7) 


R= - (19) 
4X, -I EG: -yy 
其 中 ， 
T 
yw (20) 
a TAa 
yw (21) 


AP 2 是 观测 值 ;7 是 预测 值 , -1<R<1。R 值 越 
趋 近 于 1, 则 说 明 预 测 值 与 观测 值 之 间 拟 合 越 好 。 
1.4.2 Nash-Sutcliffe 系数 (NSE) 


30r 
Xiin-»» 
NSE 一 般 用 以 验证 水 文 模型 模拟 结果 的 好 坏 ， 
取 值 为 负 无 穷 至 1。NSE 接近 1 ,表示 模拟 质量 好 ， 
模型 可 信和 度 高 ;NSE 接近 0 ,表示 模拟 结果 接近 观测 
值 的 平均 值 水 平 , 即 总 体 结果 可 信 , 但 过 程 模拟 误差 
大 ;NSE iré] T 0 , 则 模型 是 不 可 信 的 。 
1.4.3 Ai (MSE) 


NSE -1 - (22) 


1 N TS 
MSE = y X Ori - Y) (23) 


MSE 越 接近 于 0, 表 示 两 样本 (预测 值 与 观测 
值 ) 的 接近 程度 越 高 。 
1.4.4 Kendall 秩 相 关系 数 (7) Kendal 秩 相 关系 
数 是 一 种 非 参 数 检 验方 法 ,该 方法 利用 两 组 样本 
(观测 值 与 预测 值 ) 之 间 秩 的 关系 反映 二 者 之 间 的 
相关 程度 "1。 假 设 两 组 样本 为 XY,,1<i<N, 当 i 


干 区 
zjW.e 
t co I 
PXA YoY) 40 (X,-X)(Y,-Y)-20 — (Q4) 
-1 degens 
则 ， 
i up AX Y Y) 
7- NND Q5) 
7 越 接近 1 ,说 明 预 测 值 与 观测 值 之 间 拟 合 越 好 ,二 
者 之 间 相 关 程 度 越 高 。 
1.4.5  Kolmogorov-Smirnov (K-S) 检验 本文 


中 采用 两 样本 K-S 检验 ,该 方法 基于 累计 分 布 函数 ， 
用 以 检验 两 个 样本 的 经 验 分 布 是 否 相同 ,是 统计 学 
中 一 种 非常 重要 的 非 参 数 检验 方法 。K-S 统计 量 的 
计算 方法 为 : 

D = suplF,, (x) -F, „(x)l 


n,m 


(26) 


式 中 :F A F, ,分 别 表 示 两 样本 的 经 验 分 布 函数 ; 
sup 为 上 确 界 函 数 。 在 置信 水 平 w EVE 


ntm 


D, >c(a) (27) 


nm 


则 拒绝 原 假 设 。 式 (26) (27) 中 ,m,n 分别 为 两 样 
本 的 样本 容量 。 表 27 给 出 了 常用 的 置信 水 平 a 
及 其 所 对 应 的 c(a) : 

本 文 使 用 R 语言 实现 K-S 检验 ,结果 给 出 D 统 
计量 及 其 相应 的 p 值 。D 值 是 两 个 经 验 样本 之 间 的 
绝对 最 大 距离 , 称 之 为 K-S 距离 。D 值 越 接 近 0, 说 


明 两 个 样本 来 自 同一 分 布 的 可 能 性 越 大 。 
X2 置信 水 平 w 及 其 对 应 的 c(a) 
Tab.2 Value of c(w) for each level of œ 
a 0.10 0.05 0. 025 0.01 0.005 0.001 
cla) 1,22 1.36 1.48 1.63 1.73 1.95 


2 实例 应 用 


2.1 研究 区 概况 

三 江平 原 (45°01’ ~ 48?27' 56" N,130?13' ~ 
135°05'26” 玉 ,图 1) , 即 东 北平 原 东北 部 ,位 于 中 国 
东北 角 , 西 起 小 兴安 岭 东 南端 , 东 至 乌 苏 里 江 , 北 自 
黑龙 江 畔 , 南 抵 兴 凯 湖 ,总 面积 5.13 x 10* km ,三 
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Fig.1 Study area and the location of meteorological stations 


江平 原 属 温带 湿润 . 半 湿 润 大 陆 性 季风 气候 ,全 年 日 
有 照 时 数 2 400 ~2 500 h,1 月 均 温 -21~ -18 C,7 
月 均 温 21 ~22 % 。 年 降水 量 500 ~650 mm,7596 ~ 
85% 集中 在 6 ~9 月 , 雨 热 同 季 , 适 于 农业 尤其 是 优 
质 水 稳 和 高 油 大 豆 的 生长 。 尽 管 三 江平 原水 利 工程 
设施 不 断 完善 ,但 水 早 灾害 却 有 所 加 重 , 特 别 是 自 
2000 年 以 来 ,干旱 频 发 ,由 于 降水 量 持 续 减少 ,已 经 
接近 或 者 突破 历史 极 值 ,造成 了 各 种 农作物 的 减 
[e E 
2.2 数据 资料 

本 文 所 采用 的 数据 来 源 于 中 国 气象 数据 网 (ht- 
tp:/data. cma. cn) ,选取 了 三 江平 原 7 个 国家 基准 
气象 站 点 ( 表 3) 1960—2016 年 的 逐日 降水 (单位 : 
0.1 mm) 和 气温 (单位 :% ) 等 数据 资料 ,并 对 数据 进 
行 了 严格 的 修订 和 质量 控制 ,包括 采用 其 相 邻 站 点 
同期 数据 插值 进行 错误 数据 的 订正 和 遗失 数据 补漏 
及 筛选 并 剔除 无 效 数 据 二 ,确保 数据 的 采集 时 间 连 
续 .完整 ,数据 的 实 有 率 和 正确 率 均 接近 100% , 满 
足 研究 要 求 的 精度 。 


RI 气象 站 信息 


Tab.3 Information of meteorological stations 


站 台 yg 经 度 AE WR 年 平均 降 年 平均 
编号 /°E /°N /m 水 /mm 气温 /SC 
50775 Ski — 130.30 47.33 227.90 640.00 2.80 
50788 tiu 132.02 47.23 66.40 339.50 3.60 
50873 ”佳木斯 130.35 46.83 81.20 527.00 3.00 
50877 = 129.55 446.33 100.10 555.60 3.30 
50888 宝 清 132.17 46.33 83.00 548.60 3.20 
50978 鸡西 ”130.97 45.30 280.80 535.00 3.80 
50983 R 133.97 45.75 100.20 566.20 3.50 
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2.3 结果 分 析 

2.3.1 模型 拟 合 度 评 价 ”分 别 对 各 站 点 3.6.9 12 
和 24 个 月 的 SPEI 建立 ARIMA 和 ANN 模型 ,并 将 
拟 合 值 与 实际 值 进行 比较 。 通 过 RWNSE MSE 以 及 
非 参数 方法 中 的 Kendall fI K-S 距离 5 种 指标 对 模 
型 的 拟 合 度 进行 评价 验证 ,结果 见 表 4、 表 5。 

从 时 间 尺 度 来 看 ,ARIMA 和 ANN 模型 对 SPEI 
的 拟 合 情况 有 相似 的 结果 。 对 于 3 个 月 尺度 和 6 个 
月 尺度 ,ARIMA 和 ANN 模型 拟 合 的 各 站 点 的 相关 
系数 均 在 0.6 附近 ,NSE 则 都 低 于 0.5 ,说 明 ARIMA 
和 ANN 模 型 不 能 有 效 反 映 3 个 月 和 6 个 月 尺度 


表 5 ANN 模型 拟 合 度 评价 
Tab.5 ANN model fitting evaluation 


X4 ARMA 模型 拟 合 度 评 价 


Tab.4 ARIMA model fitting evaluation 


Pus 站 名 R NSE MSE Kendallr K-S 距离 
3 Siki 0.633 0.401 0.581 0.449 0.152 
富 锦 0.648 0.420 0.564 0.467 0.131 
佳木斯 0.647 0.418 0.564 0.462 0.142 

兰 0.616 0.380 0.602 0.430 0.138 

宝 清 0.631 0.398 0.584 0.441 0.130 

鸡西 0.615 0.378 0.601 0.435 0.157 

虎 林 0.626 0.392 0.592 0.440 0.138 

6 Ski 0.761 0.579 0.408 0.576 0.103 
富 锦 — 0.791 0.626 0.364 0.602 0.082 
佳木斯 0.759 0.576 0.411 0.566 0.094 

= 0.749 0.560 0.426 0.557 0.110 

宝 清 0.775 0.600 0.392 0.583 0.074 

鸡西 0.771 0.594 0.394 0.583 0.082 

虎 林 0.786 0.617 0.371 0.589 0.089 

9 稚 岗 — 0.840 0.707 0.284 0.660 0. 066 
tH 0.867 0.752 0.242 0.696 0.062 
佳木斯 0.841 0.707 0.285 0.660 0. 064 

K25 0.855 0.731 0.262 0.685 0.083 

宝 清 0.867 0.751 0.244 0.676 0.057 

鸡西 0.852 0.725 0.266 0.670 0.057 

虎 林 0.860 0.740 0.253 0.681 0.052 

12 Ski 0.934 0.872 0.124 0.785 0.035 
man 0.925 0.856 0.141 0.768 0.047 

佳木斯 0.923 0.852 0.143 0.766 0.040 

K= | 0.932 0.870 0.126 0.792 0.052 

宝 清 — 0.908 0.824 0.175 0.755 0.040 

鸡西 0.925 0.856 0.140 0.770 0.038 

虎 林 0.912 0.831 0.164 0.755 0.042 

24 Ski 0.962 0.926 0.073 0.837 0.035 
tí 0.970 0.941 0.057 0.861 0.027 
佳木斯 0.961 0.923 0.075 0.839 0.023 

(K^: | 0.965 0.932 0.066 0.851 0.037 

4$ 0.951 0.904 0.097 0.818 0.039 

鸡西 0.965 0.932 0.066 0.851 0.030 

虎 林 0.953 0.908 0.089 0.825 0.034 


/moa 站 名 R NSE MSE Kendallr K-S 距离 
3 f&[xj 0.622 0.387 0.598 0.440 0. 168 
富 锦 0.629 0.396 0.592 0.445 0.175 
佳木斯 0.621 0.385 0.597 0.427 0. 169 
œ% 0.601 0.371 0.611 0.407 0.173 
宝 清 — 0.625 0.391 0.595 0.430 0.151 
鸡西 0.615 0.379 0.601 0.431 0. 147 
虎 林 0.617 0.380 0.605 0.418 0. 159 
6 #5 | 0.756 0.572 0.414 0.567 0. 114 
富 锦 — 0.769 0.591 0.396 0.577 0. 110 
佳木斯 0.744 0.554 0.431 0.553 0. 103 
[K^ 0.742 0.551 0.435 0.546 0.127 
宝 清 — 0.761 0.579 0.409 0.559 0. 102 
鸡西 0.740 0.548 0.436 0.550 0. 105 
虎 林 0.756 0.571 0.417 0.559 0. 107 
9 钉 岗 — 0.836 0.699 0.291 0.669 0. 081 
富 锦 — 0.863 0.744 0.249 0.687 0. 080 
佳木斯 10.839 0.703 0.285 0.655 0. 092 
依 兰 0.861 0.741 0.252 0.685 0.087 
宝 清 0.862 0.743 0.251 0.676 0. 085 
鸡西 0.853 0.723 0.262 0.671 0. 093 
虎 林 0.850 0.722 0.267 0.669 0. 066 
12 #5 — 0.928 0.862 0.133 0.785 0.059 
富 锦 — 0.922 0.851 0.146 0.766 0. 070 
佳木斯 0.919 0.844 0.151 0.762 0. 061 
ik= 0.932 0.869 0.127 0.797 0. 054 
宝 清 — 0.912 0.832 0.166 0.758 0. 066 
鸡西 0.920 0.846 0.147 0.771 0.047 
虎 林 0.907 0.822 0.171 0.754 0. 052 
24 #5 — 0.963 0.927 0.072 0.840 0. 043 
富 锦 — 0.972 0.944 0.055 0.865 0. 037 
佳木斯 0.961 0.924 0.074 0.840 0. 036 
ik= 0.965 0.931 0.067 0.850 0. 059 
宝 清 — 0.952 0.906 0.094 0.814 0. 050 
鸡西 0.964 0.929 0.069 0.851 0. 043 
虎 林 — 0.953 0.908 0.088 0.826 0. 044 


SPEI 序列 的 波动 变化 ,适用 性 较 差 。 而 对 于 9.12 
和 24 个 月 的 SPET, 特 别 是 对 SPEI -12 fü SPEI -24 
的 预测 ,R 均 达 到 0.9 ME, NSE 达 0.8 以 上 ; 非 参 
数 指标 中 的 Kendallz 最 高 为 0.865 ,K-S 距离 均 小 于 
0.1, 说 明 ARIMA 和 ANN 模型 在 三 江平 原 长 时 间 尺 
度 的 干旱 中 适用 性 更 强 。 图 2 分 别 给 出 ARIMA 和 
ANN 模型 对 7 个 站 点 SPEI 拟 合 结果 各 评价 指标 的 
平均 值 。 发 现 , 随 着 时 间 尺 度 的 增加 , R, NSE Ken- 
dallz 这 三 个 衡量 预测 值 与 真实 值 一 致 性 的 指标 递 
增 ,而 MSE K-S 距离 这 两 个 衡量 预测 值 与 真实 值 不 
一 致 性 的 指标 逐渐 减 小 ,说 明 SPEI 拟 合 值 与 实际 值 
越 来 越 接 近 ,模型 的 拟 合 度 逐 渐 提 高 。 
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图 2 ARIMA 和 ANN 模型 对 7 个 站 点 拟 合 度 评价 的 平均 值 
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从 两 模型 适用 性 角度 来 看 ,对 于 同一 时 间 尺 度 ， 
ARIMA 模型 的 拟 合 效果 比 ANN 模型 略 优 ( 图 3)， 
但 这 种 优势 随 着 时 间 尺 度 的 增加 而 变 的 更 弱 。 即 整 
体 上 二 者 对 于 三 江平 原 较 长 时 间 尺 度 (9、12 .24 个 
HOSPEI 的 拟 合 效果 都 比较 好 。 例 如 对 于 12 个 月 
尺度 ,ARIMA 对 各 站 点 SPEI 拟 合 度 评价 指标 的 平 


MSE Kendallr K-S 距离 


[&] ARIMA [ANN 


图 3 不 同时 间 尺 度 两 模型 的 拟 合 度 比较 
Fig.3 Comparison of fitting for two models in different 


time scales 


均值 显示 ,RNSE MSE 以 及 非 参 数 方法 中 的 Ken- 
dallr 和 K-S 距离 分 别 为 0. 923,0. 852,0. 145, 
0.770 .0.042 ,ANN 模型 的 对 应 取 值 分 别 为 0. 920、 
0. 847 .0. 149 .0.770 .0.058 。 

2.3.2 模型 预测 精度 分 析 ”利用 最 终 确定 的 ARI- 
MA 和 ANN 模型 对 三 江平 原 7 个 站 点 2016 年 各 尺 


BARE: ARIMA 和 ANN 模型 的 干旱 预测 适 


度 SPEI 进行 预测 ,将 预测 值 与 实际 值 通过 精度 1 - 
nc 


" 行 比较 ,结果 见 表 6 ~8。 表 6 ~8 
中 仅 列 出 精度 在 70% 以 上 的 结果 。 对 于 3 个 月 和 6 
PH NEE, ARIMA 和 ANN 模型 的 预测 精度 均 在 


6096 以 下 ,表明 在 三 江平 原 地 区 两 模型 均 不 适用 于 
短 时间 尺 度 的 干旱 预测 。 对 于 9、12 .24 个 月 尺度 的 


x100% 进 


表 6 9 个 月 尺度 SPEI 预测 精度 / 96 


舌 用 性 研究 


SPET, 随 着 时 间 尺 度 的 增加 ,两 模型 的 预测 能 力 逐 
渐 提 高 。 特 别 是 对 于 12,24 个 月 SPEI 的 1 步 预 测 ， 
精度 达到 80% 以 上 。 从 模型 角度 来 看 ,对 于 所 有 的 
时 间 尺 度 , 相 比 于 ARIMA 模型 ,ANN 模型 预测 精度 
更 高 。 


R7 12 个 月 尺度 SPEI 预测 精度 / % 
Tab.7 Forecasting precision for SPEI -12 / % 


候 干旱 指数 SPEL, RI ARIMA 和 ANN 模型 对 三 江 
平原 不 同时 间 尺 度 的 SPEI 序列 建 模拟 合 预测 ,并 集 
合 参 数 和 非 参 数 方法 中 的 5 个 评价 指标 进行 对 比 验 
证 ,最 后 进行 了 12 步 预 测 ,得 到 如 下 结论 : 
从 时 间 尺 度 看 ,ARIMA 模型 和 ANN 模型 对 

SPEI 的 预测 能 力 都 随时 间 尺 度 的 增加 而 逐渐 提高 。 
对 于 3 个 月 和 6 个 月 尺度 ,预测 精度 较 低 ;对 9、12、 

24 个 月 的 SPEI 值 ,两 模型 均 可 进行 精度 在 70% 以 
上 的 1.3.4 步 预测 。 即 整体 上 ARMA 模型 和 ANN 
模型 在 三 江平 原 较 长 时 间 尺 度 (9 12 .24 ^H) FE 
预测 中 都 显示 出 良好 的 效果 。 根 据 SPEI 的 计算 公 
式 可 知 ,长 时 间 尺 度 的 SPEI 集合 了 原始 数据 中 更 多 
的 信息 ,使 得 预测 值 对 观测 值 的 拟 合 更 充分 。 


尺度 站 名 ARIMA ARIMA ARIMA ANN ANN ANN 
p 
Tab.6 Forecasting precision for SPEI -9 / 96 / mon 1] 2H 3H 1 月 2 月 34 
KE / mon "y" ARIMA 1 H ANNI 12 fp 87.36 80.56 74.75 87.73 87.26 86.11 
9 m 79.18 89.29 富 锦 90.93 82.88 73.22 $88.57 97.31 95.66 
富 锦 76. 82 94.01 佳木斯 86.89 84.06 80.48 87.35 84.98 80.70 
佳木斯 77.90 92. 64 K% 86.44 79.46 72.32 86.61 88.43 85.94 
农 兰 81.35 90. 36 E 
48 87.53 83.23 73.64 88.78 85.01 81.55 
EX. 82.21 82.25 
鸡西 84. 50 89. 56 35 88.66 77.28 77.92 090.00 81.77 86.07 
虎 林 72. 19 81.93 虎 林 83.71 76.45 70.93 85.04 86.03 86.75 
表 8 24 个 月 尺度 SPEI 预测 精度 / 96 
Tab.8 Forecasting precision for SPEI —24 / 9o 
尺度 站 名 ARIMA ARIMA ARIMA ARIMA ANN ANN ANN ANN 
/ mon a 1H 2 月 3 月 4 H 1 月 2 月 3 月 4 H 
24 稚 岗 90.23 85.09 83.95 80. 19 93.47 89.34 87.69 83.80 
富 锦 95.07 85.32 96.58 78.44 95.36 92.15 95.63 86.12 
佳木斯 85.70 87.83 82.56 76. 97 92. 98 89. 84 87.97 83.02 
依 兰 96.21 92.86 91.80 92. 84 95.39 92.23 96.28 96.36 
a 88.00 84.54 71.66 70.98 89.13 85.35 73.05 70.03 
鸡西 90.72 85.22 82.37 82.50 94.75 87.74 89.23 86.39 
虎 林 82.69 74.79 75.71 75. 43 85.20 83.71 83.96 77.70 
4t 论 对 于 相同 时 间 尺 度 , 在 训练 和 测试 阶段 , ARI- 
£n VE 
MA 模型 的 拟 合 效果 与 ANN 模型 基本 持平 ,但 是 在 
本 文 引 入 了 综合 降水 和 气温 变化 共同 效应 的 气 。” ”预测 时 ,ANN 模型 则 表现 出 相对 更 高 的 精度 。 说 明 


基于 ANN 时 间 序 列 预测 模型 具有 较 强 的 可 操作 性 
可 以 有 效 地 进行 三 江平 原 较 长 时 间 尺 度 干旱 情况 的 
预测 。 这 主要 是 由 于 以 下 两 方面 原因 导致 :(1) 降 
雨 数据 的 产生 是 随机 的 .不 确定 的 ,并 且 大 多 数 数据 
都 含有 白 噪 声 , 而 ARIMA 模型 本 质 上 是 一 种 整体 
线性 自 回 归 模 型 中 ,对 于 由 降水 数据 计算 得 来 的 
SPEI 的 变化 这 种 非 线性 行为 的 分 析 、 预 测 存在 着 较 
大 的 不 确定 性 和 不 足 ;ANN 则 具有 和 较 强 的 非 线 性 映 
射 能 力 和 较 好 的 稳健 性 '" ,能 够 区 分 有 噪声 的 样 
本 ,可 以 通过 对 样本 学 习 把 隐 含 的 内 在 规律 分 布 在 
网 络 的 连接 权重 上 。(2) ARIMA 模型 的 基本 思想 
是 首先 将 非 平 稳 的 时 间 序 列 通过 差分 转化 为 平稳 的 
时 间 序 列 ,再 用 ARMA 模型 对 该 平稳 序列 建 模 ”， 
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而 在 差分 的 过 程 中 不 可 避免 的 会 损失 数值 信息 ; 
ANN 模型 则 充分 发 挥 了 较 强 的 学 习 能 力 和 通 近 任 
意 函 数 的 特点 ,对 原始 数据 进行 了 充分 的 学 习 。 

此 外 ,对 于 短 时 间 尺 度 的 干旱 预测 有 待 在 日 后 
的 研究 中 完善 。 
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Applicability of ARIMA and ANN models for drought forecasting 


YANG Hui-rong', ZHANG Yu-hu’, CUI Heng-jian!, GAO Feng’, CHEN Qiu-hua' 
(1 School of Mathematical Sciences , Capital Normal University , Beijing 100048 , China ; 
2 College of Resources Environment & Tourism , Capital Normal University , Beijing 100048 , China.) 


Abstract: Drought is one of the major natural disasters , whose occurrence is linked to a sustained lack of precipi- 
tation. The drought forecast provides vital evidence and support for preventing losses of drought disasters , and there- 
fore it is of great significance. In this study ,a series of the standard precipitation evapotranspiration index ( SPEI) at 
different time scales were calculated based on the daily precipitation and temperature data from 7 meteorological sta- 
tions in Sanjiang Plain , northeast China from 1960 to 2016 and were used to forecast the drought using ARIMA and 
ANN models. In the stage of training and testing ,the fitting degrees of the models were evaluated and validated and 
the optimal ARIMA and ANN models were chosen with the help of 6 fitting evaluation methods; the correlation coef- 
ficient ( R) , Nash-Suicliffe efficiency coefficient ( NSE) , Kendall , rank correlation coefficient ,the mean square error 
( MSE) and Kolmogorov-Smirnov ( K-S) test. Then 12 values for the 12 months in 2016 were predicted by the opti- 
mal models and were compared with the corresponding observations. The results are shown as follows: (1) The pre- 
diction ability of ARIMA and ANN models based on SPEI were both increased with the increase of time scale in 
Sanjiang Plain. (2) The two models had poor prediction accuracy for SPEI 3 and SPEI 6. For the SPEI value of 9, 
12 and 24 months , all models worked well with accuracy more than 7096. (3) For the SPEI value of 9,12 and 24 
months , the prediction accuracy of ANN model is better than that of ARIMA model. In particular,the prediction ac- 
curacy for one month forecast of SPEI 12 and 24 at all stations were more than 80% . All these showed that the pre- 
diction model of ANN has strong maneuverability and can effectively predict the drought at a large time scale in San- 
jiang Plain. The drought prediction at small time scale (3 and 6 months) needs to be improved in future studies. 


Key words: drought; ARIMA model; ANN model; SPEIs; Sanjiang Plain 


